Spin-orbit induced anisotropy in the magneto-conductance of two-dimensional metals 
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It is shown that the spin-orbit coupling due to structure inversion asymmetry leads to a char- 
acteristic anisotropy in the magnetoconductance of two-dimensional metals. Relevance for recent 
experiments is discussed. 
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Over the last five years there has been a signifi- 
cant amount of experimental evidence of a metallic-like 
temperature dependent resisthdty for a certain electron 
density range in Si- MOSFETsU and semiconductor het- 
erostructures. This has called for a better understanding 
of the possible ground states of two-dimensional (2D) 
systems. For such an endeavour it is important from 
the outset to identify the relevant physical mechanisms 
at play. A recent summary of the main experimental 
facts, and suggested theoretical models can be found in 
Ref.EI. In the 2D systems which are experimentally inves- 
tigated, the spin-orbit interaction, as due to the lack of 
structural inversion symraetryu, has been proposed as a 
potentially relevant mechanism for the observed resistiv- 
ity behaviouro. However, there is no consensus at present 
about the importance of this spin-orbit interaction. 

Magneto-transport experiments are very powerful 
tools in selecting processes which are affected by the pres- 
ence of a magnetic field. In the case of 2D systems one 
of the most puzzling features is the magnetoresistance in 
a parallel magnetic field. For carrier densities of interest 
the resistance increases by roughly one order of magni- 
tude on the scale of several Tesla before saturating at 
high magnetic fields! It has been generally assumed that 
a parallel magnetic field affects the system only via the 
coupling to the electron spin, and there is experimental 
evidence supporting this view. However the importance 
of the coupling of the magnetic field to the orbital mo- 
tion of the electrons as a consequence of the finite layer 
thickness of the 2D structure has also been proposed!! 

In order to shed light on the problem there has re- 
cently been the proposal that, in the presence of the spin- 
orbit interaction, an in-plane magnetic field gives rise to 
an anisotropy in the magnetoresistance, depending on 
whether one measures the current in the direction per- 
pendicular or parallel to the magnetic field! This has 
stimulated more experimental work, where anisotropy 
in the magnetoresistance has been found in GaAs hole 
systems!, Si-MOSFETsfJ, and GaAs electron systemsM 
The prediction of Ref.Q has been made in the deeply insu- 



lating limit and this makes a direct comparison with the 
experiments a little difficult. This fact has motivated us 
to investigate the interplay of the spin-orbit interaction 
and an in-plane magnetic field in the metallic regime. 
Here we calculate the conductivity in the framework of 
the Drudc-Boltzmann theory. We do not worry about 
the microscopic orign of the scattering which is, for ex- 
ample, responsible for the strong magnetic field depen- 
dence of the conductivity. Our concern is the anisotropy 
of the conductance and its behaviour as a function of the 
various physical parameters. The most important find- 
ings are: a) the magnetoresistance is anisotropic; b) the 
sign of the anisotropy depends on the type of scattering 
potential seen by the charge carriers; c) the anisotropy 
factor has a maximum as function of the magnetic field, 
with the maximum position scaling with the density. In 
the following we present details of our calculations and 
discuss its possible experimental relevance. 
We start from the model Hamiltonian 
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with a a parameter describing the spin-orbit interaction 
due to the confinement field and a being a vector whose 
components are the Pauli matrices. The unit vector e 2 
is perpendicular to the 2D plane and defines the z-axis. 
The magnetic field is chosen to lie in x-direction. One 
finds then two bands with the dispersion 



E±(p) =p 2 /2m±n(p), 



(2) 



with Q(p) = |a:p A e 2 — ui s e x \. 

In the problem we have several energy scales: the 
Fermi energy tp (in the absence of magnetic field with 
no spin-orbit coupling), the Zeeman energy ui s = ^gusB 
(note that this differs by a factor two from the stan- 
dard definition), and the spin-orbit energy apF- It is 
useful to make a few estimates. In Si inversion layers one 
has a density-of-states Ao = 1.59 x 10 11 cm _2 mgV _ , 
which for the range of densities considered in Reffl (i.e. 
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n = 0.7 — 1.3 x 10 n cm~ 2 ) gives for the Fermi energy 
ep — 5.1 — 9.5K. By assuming g — 2 a magnetic field of 
1 T gives uj s =p 0.6K. For the spin-orbit interaction, we 
take from RefB a — 6 x 10~ 6 K cm, which gives a spin- 
orbit energy apF — 3.36 K. Considering transport intro- 
duces another scale, which we characterize by the scat- 
tering rate 1/r. We assume throughout this paper that 
epT is larger than unity as is appropriate for the metal- 
lic regime (in units with % = 1). A realistic mobility of 
4 x 10 4 cm 2 /Vs gives then a relaxation rate 1/r <~ 1 — 2K. 
In the following we choose parameters having the Si sys- 
tems in mindtll. Measuring energies in units of €f one has 
that 1/t = 0.1-0.4, ap F = 0.3-0.7 and u s = 0.06-0.15 
per Tesla. 

A simple estimate for the conductivity and its 
anisotropy may be obtained with the semi-classical Boltz- 
mann equation in the relaxation time approximation. As- 
suming a single transport time, i.e. momentum- and 
band-independent, the conductivity is 
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where the sum is over the two bands, N± is the density 
of states, v± = V p E±(p) the velocity, and (. . .)fs the 
average over the Fermi surface. It is apparent that the 
transport time drops, when we consider the conductivity 
relative to its value in the absence of magnetic field and 
spin-orbit coupling, ctq = e 2 Nov F Tt r - 
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FIG. 1. Anisotropy in the magneto-conductance as ob- 
tained from the Boltzmann equation in the relaxation time 
approximation. 

In Fig.|l| we plot the anisotropy in the magneto- 
conductance for various strengths of the spin-orbit cou- 
pling. One observes a negative anisotropy, i.e. the con- 
ductivity perpendicular to the field is larger than that 
parallel to the field, with a maximum in the region 
lu s ~ €f- Some insight in the results reported in Fig.|l| 
is obtained from the analytic expressions in the weak 
(u> s < ep) and strong magnetic field limits. For the given 
dispersion relation the velocities are v± = p/m ± ae p , 
with e p = (ap — w s e y )/Q(p). In the weak magnetic field 
limit, one finds that only the anomalous part of the ve- 
locity contributes to the anisotropy in the conductance. 
We obtain then 

<Jxx - o-yy = 2e 2 N T tI a 2 ((e p ■ e x ) 2 - (e p • e y ) 2 ) FS . (4) 



For apF > u s the polarization vector e p makes a full 
rotation when averaging over the Fermi surface. In the 
limit lo s > apF however the polarization vector becomes 
locked with e p w -e 9 . It is clear that one then finds a 
negative anisotropy in the conductance. Going explicity 
through the algebra we arrive at 
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which has an edge at uj s — ap F and becomes constant 
when uj s ^> apF- In the strong field limit the arguments 
are different. For large field only the lower band is occu- 
pied. The Fermi surface again is almost rotational sym- 
metric, but there is a weak elliptic distortion which leads 
to an anisotropy: 
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From these considerations we find that the anisotropy for 
ui s r~j €f scales as (apF Af) 2 ■ 

In order to have a more microscopic calculation of the 
conductivity we now go beyond the relaxation time ap- 
proximation. We introduce a disorder potential U, which 
will be treated in the self-consistent Born approximation 
and the conductivity will be calculated using the Green 
function formalism. For simplicity we assume a impurity 
potential with short range Gaussian correlations 



( C /(x)C/(x')) = -4 r -5(x-xO. 
zttNqt 



(7) 



The self-energy is determined self-consistently by the 
equation 
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where G(p) is the Green function. Notice that S and 
G(p) are 2x2 matrices. Expanding the self-energy in 
Pauli matrix components, it turns out that the self- 
energy has no 02 and 03 components. As a result the 
self-energy will have the form S = Eoco + ^i&i- The 
real part of the self energy So shifts the energy spectrum 
by a constant. Since we have to adjust the chemical po- 
tential in order to keep the particle number at a given 
value Re£ can be safely neglected. The real part of £1 
gives rise to a renormalization of the Zeeman energy. For 
the imaginary part, we find in the limit of weak disorder 

Im£^ = -l/2r « -(N+ + N_)/4N q t (9) 
Im£f = -1/2ti w -(N- - A+)/4A r (10) 

The sum or difference l/r T = l/r ± 1/n are roughly the 
scattering rates for the two subbands. For weak magnetic 
field (lu s < ep), the density of states in the two subbands 
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are identical and therefore l/r± = 1/to = 1/r. In the 
strong field limit (w s > tp) the upper band is depopu- 
lated, so that 1/ti = 1/t = l/2r. 

The current operator, which is necessary for the eval- 
uation of the conductivity, is modified in the presence of 
the spin-orbit interaction: 
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The conductivity is obtained by the formula 

= £ E Tr [fG R JLG A - jW^G 1 
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The dressed vertex J 3 depends upon whether it is con- 
nected to a pair of retarded and advanced Green func- 
tions or a pair of Green functions with equal analytic 
properties. 

In order to evaluate the conductivity, we have to deter- 
mine the renormalized vertices and perform the momen- 
tum integrals. Before showing the complete results it is 
useful to consider the conductivity in the absence of ver- 
tex corrections and compare with the above Boltzmann 
equation analysis. 
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FIG. 2. Anisotropy in the magneto-conductance for 
app = OAep, l/(€pr) = 0.4; the full line is the result from 
the Boltzmann equation. The dashed lines result from the 
Kubo formula, when vertex corrections are neglected, but 
using the self-consistently determined self-energies in the 
Green functions. 

The result is shown in Fig. (|^) (dashed lines). The con- 
ductivity may be expressed as a sum of intra- and in- 
terband contributions. Because of the gap between the 
two spin subbands, the interband terms do not contribute 
much unless the disorder is strong enough to produce a 
broadening of the Fermi surface of the two bands. The 
condition for weak disorder is then Sir ^> 1. In the weak 
disorder limit one reproduces - when lu s < ep — the re- 
sults obtained from the Boltzmann equation (the solid 
line). For large magnetic field the Boltzmann result is 
not reproduced, even for very weak disorder. The reason 
is that even an isotropic scattering potential generates an 
anisotropic quasi-particle lifetime as demonstrated here 
below. Let us denote the scattering probability from 
state p to p' by W PP ' . When going from the spin- to the 



eigenstate basis the scattering probability in the lower 
band becomes 



W pp , -> FVpp'i (1 + e p • e p -) « W pp , [ 1 
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From this anisotropic scattering probability the inverse 
lifetime of an electron in state p is determined as 1/t p w 
[1 — (ap x ) 2 / 4uj^} / t which definitely depends on the posi- 
tion on the Fermi surface. The anisotropic scattering rate 
then also contributes to the anisotropy in the magneto- 
conductance. 

We are now ready to consider the vertex corrections. 
For the choice of disorder potential we have made, these 
corrections usually vanish expressing the fact that the 
relaxation of momentum is governed by the quasiparti- 
cle lifetime. However, in the present problem there are 
specific current vertices due to the spin-orbit interaction 
which acquire a vertex correction. The vertex corrections 
are obtained in the standard way by solving the equation 
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By separating the p-dependent (i.e., proportional to the 
p vector) and p- independent parts f = p 1 /m + Y and 
J 1 = Jg(p) + r l , one notices that as a consequence of 
isotropic scattering, the p-dependent part is not dressed, 
i.e., J 1 = p l /m + r\ The equation for the momentum 
independent part of the current vertex reads 
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We included here the matrix (=spin) indices s, s', a, and 
b. The quantities 7* are a sum of the bare vertices 7* and 
a term which is generated by p/m, 
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The above equation may be solved by expanding all ma- 
trices in terms of the Pauli matrix basis, e.g., T l — 
J2p,=o 3 °>rji- By skipping the details of the derivation 
we present the numerical results for the renormalized ver- 
tices and the conductivity in Figs.@ and (||), respec- 
tively. One observes that for weak magnetic field the 
quantities T become very small so that Jra ~ p/m i.e. 
the anomalous velocity is cancelled by the vertex correc- 
tions! Only for fields which are at least of the order of 
the Fermi energy does an anomalous velocity contribu- 
tion survive. 

Concerning the conductivity the most striking result is 
that the vertex corrections, for not too strong magnetic 
field, change the sign of the anisotropy. (See Fig.0). 
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FIG. 3. The dressed vertices Fq . . . T3 in units of a as a 
function of the magnetic field. Here we chose apF = OAtF 
and 1/(£ft) = 0.4. Remember that the bare vertices are 
7o = 73 = and 71 = —72 = a; the asymptotic values in 
the strong field limit are Ti = —To = a/2, T2 = —a, and 
F3 = 0. Notice that T2 reaches the asymptotic value for high 
magnetic field only very slowly. 
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FIG. 4. Anisotropy in the magneto-conductance for 
l/(epr) = 0.2 and various spin-orbit energies. For com- 
parison we included the results of the Boltzmann theory 
(apF = 0.4ef; full line). For not too strong fields, the vertex 
corrections change the sign of the effect. 

To understand the effect of the vertex correction, it is 
again useful to consider the scattering rate in the eigen- 
state basis, see eq.(|l^). One then observes that along the 
x direction, forward scattering is enhanced compared to 
backward scattering, while along the y direction there is 
no difference between forward and backward. Because 
vertex corrections increase conductivity if forward scat- 
tering is favoured, one expects that vertex corrections en- 
hance a xx most contributing to Act with a positive term, 
which for low and moderate magnetic fields overcomes 
the negative contribution discussed previously. 

We now comment on our results in the light of the 
available experiments. Both in the relaxation time ap- 
proximation and in the Green function calculation the 
maximum anisotropy in the conductivity occurs at uj s ~ 
6f irrespective of the value of the spin-orbit energy a. 
This means that the energy of the peak scales with 
the electron-, density. This is in agreement with the 
experimental A second feature is that the peak strength 
is controlled by the spin-orbit energy in units of ep. 
Hence the peak strength must scale with the inverse den- 
sity, so that the effect is expected to decrease going more 
deeply into the metallic regime. Thia_again is in agree- 
ment with the experimental findingSLlUj. Concerning the 
sign of the effect, inclusion of the vertex corrections in 
the Green function calculation spoils the agreement with 



the experiments. A possible reason for this discrepancy 
may be due to our simplifying assumption of a short 
range p-independent disorder potential. Our analysis has 
shown in fact that the sign of the anisotropy depends 
crucially on the effective p-dependence of the scattering 
rate in the eigenstate basis. It is then reasonable to ex- 
pect that the specific choice of the scattering potential 
may indeed play a role in determining the sign of the 
effect a low fields. This analysis, while it is worth pur- 
suing, is however beyond the scope of the present pa- 
per. The anisotropy in the magnetoresistance observed 
in GaAs heterostructuresH presents a rather complex pat- 
tern. Two crystallographic directions with different mo- 
bilities have been investigated. The one with higher mo- 
bility, shows anisotropy with both Act < at not too 
large magnetic fields and Act > at larger fields. At low 
fields (where Act < 0) the value of the field at which the 
anisotropy is maximum seems to scale with the density 
(cf. right column of Fig. 2 of Ref.Q). At this stage one 
can only speculate that the actual sign of the anisotropy 
may depend on the interplay of the details of the scatter- 
ing potential with the spin-orbit coupling effect discussed 
here. 

In conclusion, our theory explains how the spin-orbit 
coupling may give rise to the anisotropy in the conductiv- 
ity, although the interpretation of the experiments may 
require more realistic models of disorder, as for example 
that due to scattering by impurities located outside the 
2D plane, which favours small-angle scattering. It is also 
our hope that this work will stimulate more experimental 
effort toward the investigation of anisotropic conductiv- 
ity as a function of the scattering potential or in sample 
with different mobilities. 
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